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We discuss lepton flavor violation (LFV) associated with tau leptons in the general framework of the two Higgs 
doublet model, in which LFV couplings are introduced as a deviation from Model II Yukawa interaction. Parameters 
of the model are constrained from experimental results and also from requirements of theoretical consistencies such 
as vacuum stability and perturbative unitarity. Current data for LFV rare tau decays provide substantial upper 
limits on the LFV Yukawa couplings in the large tan /3 region, which are comparable with predictions in fundamental 
theories. Here tan is the ratio of vacuum expectation values of the two Higgs doublets. A search for the LFV decays 
0° — > r^^T (7±gT) of neutral Higgs bosons (<f>° = h,H and A) at future collider experiments can be useful to further 
constrain the LFV couplings, especially in the relatively small tan j3 region (tan j3 < 30), where rare tau decay data 
cannot give any strong limit. 



1. Introduction 



> : 

The structure of the electroweak symmetry beaking sector would directly connect to the property of physics beyond 
the standard model (SM). Models of such new physics predict extended Higgs sectors with more than one scalar 
doublets in the low energy effective theory. These extended Higgs models would show distinctive features from the 
, SM phenomenology. The most obvious evidence would be the confirmation of the existence of the extra scalar states 
such as CP-odd and charged states at future collider experiments. Even if they are too heavy to be directly discovered 
and only the lightest Higgs boson is found at the experiments, we would be able to explore the extended Higgs sector 
by searching for deviations from the SM predictions on the couplings with gauge bosons and fermions as well as on 
the self coupling. Furthermore, search for non-SM interactions can also be useful. 

Lepton flavor violation (LFV) is an example for such non-SM phenomena. In particular, LFV in the Yukawa sector 
can only appear for extended Higgs sectors. Flavor violation between electrons and muons 1] has been tested through 
rare muon decays such as \i — ► ej and fi — > ee+e~, as well as through u-e conversion. Tau lepton associated LFV 
has also been studied by rare decays of tau leptons such as r -> ^P°H) t ^M+M' - 0,3) T ~> 
and r — » ^70,0,0, where li (i — 1, 2) respectively represent an electron and a muon, P° does tt° , r] and rf mesons, 
(M ,=1= ) does 7r ± and mesons, and £ ,ziz = e ± and fi ± . The LFV Yukawa couplings can be constrained from 
the data for these processes especially those with the Higgs boson mediation. For /i-e mixing, the Higgs boson 
mediated LFV coupling has been discussed in Ref . [lfj, . Tau lepton associated LFV processes with the Higgs 
boson mediation have been discussed in models with s upersym metry fSUSY)[l2l \l$L 0] as well as in the two 
Higgs doublet model (THDM) in some specific scenarios 

liaLUllia- In Ref - 113 > tau associated LFV processes have 
been discussed comprehensively in the framework of 4- Fermi contact interactions. Phenomenological consequences 
of the LFV Yukawa couplings associated with tau leptons have also been studied for future observables such as B s 
decays [HQ at (super) B factories [il] and Higgs boson decays El El at CERN LHC0, an electron-positron 
linear collider (LC) 26] and a muon collider |27j. Furthermore, deep inelastic scattering processes fxN — > rv|2^. Eflf 
from intense high energy muons at neutrino factories (or muon colliders) and eN — ► tY 29] by using the electron 
(positron) beam of a LC would be useful to further explore the tau lepton associated LFV Yukawa couplings. 
In this talk, we discuss LFV in Higgs boson decays into a T-li pair in the general framework of the THDM0]. We 
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study experimental upper limits on the tau lepton associated LFV Yukawa couplings to evaluate possible maximal 
values of the branching fractions. The parameter space of the model is tested by theoretical requirements for vacuum 
st ability and perturbative unitaritv ^llE^ . Current data from electroweak precision measurements at LEpj^Ul^j 
and those at the B factories |35| also strongly constrain parameters of the Higgs potential. Under these theoretical 
bounds and experimental limits on the model, possible maximal values of the LFV couplings of r-^j-0 are obtained 
by using the current data for rare tau decays, where 0° represents two CP-even (h and H) and a CP-odd (A) Higgs 
bosons. We then evaluate branching ratios of 0° — ► r ± £? with the maximal allowed values of the LFV couplings of 
T-£i-(j)° in a wide range of the parameter space. 



(1) 



2. Lepton flavor violation in Yukawa interaction 

The Higgs sector of the general THDM is expressed as 

-£ffi ggs = m 2 |$i| 2 +m 2 |$ 2 | 2 - (m^^+H.c.) +^|$ 1 | 4 + y |$ 2 | 4 + A 3 |$i| 2 |$2| 2 + A 4 |$ t 1 $ 2 ' 2 

+ |y ($ t 1 $ 2 ) 2 + H.c.| + (\ 6 |$i| 2 $J$ 2 + H.c.) + (a 7 |$ 2 | 2 $i$2+H.c 

where $i and $2 are the scalar iso-doublets with hypercharge 1/2. In Eq. to 2 , A5, Ag and A7 are complex in 
general. We here assume that all the parameters to 2 _ 3 and Ai_7 are real. The terms of to 2 , A6 and A7 break the 
discrete symmetry explicitly. As we consider the model in which the discrete symmetry is explicitly broken only in 
the leptonic Yukawa interaction, we set the hard-breaking coupling constants to be zero in the Higgs potential; i.e., 
A 6 = A7 = 1 , and retain only the soft-breaking mass parameter m§. 

There are eight degrees of freedom in the two Higgs doublet fields. Three of them are absorbed by the weak 
gauge bosons via the Higgs mechanism, and remaining five are physical states. After the diagonalization of the mass 
matrices, they correspond to two CP-even (h and H), a CP-odd (A), and a pair of charged (H^) Higgs bosons. We 
define such that h is lighter than H. The eight real parameters to 2 _ 3 and Ai_5 can be described by the same number 
of physical parameters; i.e., the vacuum expectation value v (~ 246 GeV), the Higgs boson masses m h , m H , m A and 
m H± , the mixing angle a between the CP-even Higgs bosons, the ratio tan/3 (= ($ 2 ) / (^1)) of the vacuum expectation 
values for two Higgs doublets, and the soft-breaking scale M (= y/ m\ / sin (3 cos /3) for the discrete symmetry. The 
quartic couplings are expressed in terms of physical parameters in Ref. 

Parameters of the Higgs sector are constrained from requirements of theoretical consistencies and also from the 
current experimental results. We here take into account two kinds of theoretical conditions; i.e., vacuum stabilitv[3oj 
and perturbative unitaritv|3ll 1 3 2j | at the tree level. The condition for tree-level unitarity, which we employ in 
Refs. [32l l.36| . is described as | { 03 04 1 a. 1 </>i 02 ) | < where (0 3 04|a°|0i02) is the s-wave amplitude for the process of 
0102 — > 0304 with a (a—l-A) denoting Higgs bosons and longitudinal components of weak gauge bosons. We take the 
criterion £ to be 1 (and also 1/2 for comparison). The experimental constraints are prov ided by the LEP precision 
data[3|, the b — ► sj results ji^, and the direct search results for the Higgs bosons j23|. The LEP precision data 
provide the strong constraint on the new physics structure via the gauge-boson two-point functions. The constraint 
on p parameter indicates that the Higgs sector is approximately custodial SU(2) symmetric. This requirement is 
satisfied when (i) rn H± ~ rn A , (ii) m H ± ~ m H with sin 2 (a — (3) ~ 1, and (hi) m H ± ~ m h with cos 2 (a — /?) ~ 1. It 
is known that in Model II, the b — > 37 result gives the lower bound on the charged Higgs boson mass. We here take 
into account this bound by requiring m H± > 350 GeV. 

Next, we consider the Yukawa interaction for charged leptons as 

- Acpton = 1m { Yu % *i + (Y h t% + , U Y, )$ 2 }-L j + H.c. , (2) 



1 Even in such a case, Ag and A7 are effectively induced at the loop level. They are suppressed by the loop factor, so that we here 
neglect these small effects. 
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where t R i (£=1-3) are right-handed charged leptons, and Li (i=l-3) denote the lepton doublets and Yi t are the 
Yukawa couplings for l^. This interaction is reduced to be of Model II 37] in the limit ef^ R — > with the discrete 
symmetry under e R — > +e R , Li — > +Li, $i — > +<f>i, and $2 — > —$2- Nonzero values of t^ R (i ^ j) yield the LFV 
Yukawa couplings after the diagonalization of the mass matrix. We note that in supersymmetric standard models, 

L R 

the Yukawa interaction for leptons is of Model II at the tree level, and e-' can be induced at the loop level due 
to slepton mixing (13. fisl . For the quark sector, Model II Yukawa interactions are assumed to suppress FCNC, 
imposing the invariance under the transformation of u R — > — u R , d l R — ► +d R , q l L — > $1 — > +$1, and $2 — > — $2- 
The tau lepton associated LFV interactions in Eq. J5J can be reduced in the mass eigenbasis of each field to 

£ tL fv = m \ a (KaifPjJi + K%!iP L T) {cos {a-@)h + sin (a — f3) H — iA} 



w cos 2 /3 
•\/2rn T 



-r 5-5 V 4^ P i^ + <3^P^r J?" + H.C., (3) 

v cos p v ' 

where is the projection operator to the left-handed field, and l\ and i? 2 respectively represent e and /x. In general, 
the LFV parameters k^- r can be expressed in terms of e^- R and tan (3. 2 We here take these k^- R as effective couplings, 
and investigate their phenomenological consequences. We note that Eq. © is exact in the limit of m e , — > 0. The 
terms of Kf 3 and K R t (i = 1,2) are proportional to m» , so that they decouple in this approximation. 

L R 

We briefly discuss relationship between k- ■' and new physics models beyond the cut-off scale of the effective 

L R 

THDM. When a new physics model is specified at the high energy scale, can be predicted as a function of the 
model parameters. For example, in the MSSM, slepton mixing can be a source of LFV. Notice that the induced LFV 
Higgs interactions do not necessarily decouple in the limit where the SUSY particles are sufficiently heavy, because 
their couplings only depend on the ratio of the SUSY parameters. Therefore, the Higgs associated LFV processes can 
become important in a scenario with the soft-SUSY-breaking scale to susy to be much higher than the electroweak 
one. In the MSSM, predicted values of \n 3i \ 2 can be as large as of 0(1O~ 6 ) when to susy is a few Tev[liEi- In 
the MSSM with right-handed neutrinos, left-handed slepton mixing may be a consequence of running effects of the 
neutrino Yukawa couplings between the scale of the grand unification and that of the right-handed neutrinos |38f. 
The parameters K 3i are mainly induced by mixing of left-handed sleptonsfl^. lul. [Lll fljl Ejj . . The LFV Yukawa 
interactions can also appear effectively in the Zee model j^Sj. The LFV parameters are induced through flavor 
violating couplings in the charged scalar interactions with leptons. 



3. Bound on LFV Yukawa couplings from rare tau decays 

In order to constrain the LFV parameters |«^| and we take into account the data for rare tau decay processes 
such as t — > ivP , t -> liM+M'-, r -> id'+l'-, and r -> £7, where P° represents tt°, 77 and rf mesons, M ± (M f± ) 
does tt^ and K mesons, and £ /=t = and fi^. The list of the current data from the B factories are summarized 
in Table U [E IE El IE IE EL IE These bounds may be improved at the super B factory around a digit [2l|. In our 
analysis, we take the underlined data in Table [i] as our numerical inputs. The branching ratios for these rare r LFV 
decays are summarized in Ref. j^E]- 

Since branching ratios for r _ — > £^P°, t~ — > tjM + M~, t~ — > i~£' + £'~ and r~ — > £^7 depend on different 
combinations of the Higgs boson masses, independent information can be obtained for the model parameters by 
measuring each of them. When all the masses of Higgs bosons are large, these decay processes decouple by a factor 



2 LFV parameters nfj R can be expressed as 



$ = r^z (X = L,R), (4) 



{l + (ef 3 + ef 3 ) tan/3} 

for the case of e^' fl tan/3 <C O(l) which is satisfied in the case of the minimal supersymmetric standard model (MSSM~l ll2lll3l 
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Table I: Current experimental limits on branching ratios of the LFV rare tau decays. 



of l/ m Higgs- Although these branching ratios are complicated functions of the mixing angles, each of them can be 
simply expressed to be proportional to tan 6 /3 for tan/3 ^ 1 in the SM-like region (sin(a — /3) ~ —1 |3EHE3l)- 
This tan 6 f3 dependence is a common feature of the tau-associated LFV processes with the Higgs-mediated 4- Fermi 
interactions. 

The experimental upper limit on | | 2 (= Ik^I 2 + |^| 2 ) can be obtained by using the experimental results given 
in Table |U and analytic expressions of the decay branching ratios for rare tau LFV decay processes in Ref . 24] . For 
description, let us consider the bound from the r — * fir) results [l4j: 

|„. 12 / / | t .max|2\ _ 2567rBr(r ^/_ii)) cxp mj L cos 6 3 / r\ 

|«32| S I |« 3 2 I ) = 7 rjTT^rgi {^J 

\ /T^flT] Q(72_ m 3 m 4 F2 r / 1 V 



where Br(r — > nrj) cxp is the experimental upper limit on the branching ratio of r — * fj,rj in Table [I] In particular, for 
tan /?> 1, the right-hand-side can be expressed by 

(\ K ?r\ 2 ) ^2.3xl(T 4 x ( ™ A - Vf JLV. (6) 
V 132 1 Jr^u,, UsOGeV/ V tan/3/ W 



350[GeV]7 Vtan/3 

It can be easily seen that the bound (Inf^^r^^ is rapidly relaxed in the region with small tan/3 and large m A . In 
a similar way to Eq. (J3J), the maximal allowed value (l K 3i ax | 2 ) modo can be calculated for each mode. The combined 
upper limit |«^ ax | 2 is then given by 



| K £ ax | 2 = MIN 



As shown below, r — > ^ry and r — > ^7 give the strongest upper limits on \n^i\ 2 in a wide range of the parameter 
space. In addition, in some parameter regions r — > £ i K + K^ and r — > l,-^ can also give similar limits on |K3i| 2 
to those from the above two processes. 
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(a) 



tan/3 



(b) 



tan/3 



Figure 1: Contours of |kJ) 



the possible maximal value of |k32| 2 from the rare tau decay results, are shown (a) in the 



tan (i-vriA plane and (b) in the tan/3-sin 2 (a — /3) plane. The parameters are taken to be (a) m h — 120 GeV, m H — m H ± = 350 
GeV and sin(a — 0) = —0.9999, and (b) m h = 120 GeV and m H — uia = ra H ± = 350 GeV. The remaining parameter M is 
scanned from to 1,000 GeV. The dark (light) shaded area indicates the excluded region by the theoretical requirements of 
vacuum stability and perturbative unitarity with £ = 1 (£ = 1/2). 



In Figs. ID (a) and[T]-(b), contour plots for |K^ ax | 2 are shown under the rare tau decay results in the tan/3-m^ 
plane and the tan/3-sin 2 (a — 0) plane, respectively. The combined excluded region from the theoretical conditions 
of vacuum stability and pertrbative unitarity is indicated by the dark shaded area for the criterion £ = 1 and by the 



= 120 GeV, m H = m H± = 350 
= 350 GeV. 



' H 



light one for £ = 1/2. In Fig.nj(a), parameters of the Higgs sector are taken to be r 
GeV and sin(a - 0) = -0.9999. In Fig. QJ(b), those are m h = 120 GeV and 
The value of |ft§?f x | 2 is independent of M, the soft-breaking scale of the discrete symmetry. On the other hand, 
theoretical bounds from vacuum stability and perturbative unitarity are sensitive to M. Therefore, we evaluate such 
a theoretical allowed region by scanning M to be from to 1000 GeV. We also take into account the constraint from 
the p parameter measurement and the b — > S7 result by taking sin(a — /3) ~ —1 and m H = m 



H ± with m H± > 350 



GeV for Fig^-(a), and m A = m H± with m H± > 350 GeV for Fig^-(b). From the both figures, it is easily found that 
the value of |«32 ax | 2 can extensively be larger for smaller tan/3 in the allowed region under the theoretical constraints. 

? 32 ax | 2 can be O(0. 1) (O(10~ 4 )). Among the rare tau decay processes, r — > jirj and r — > hj 



For tan/3 < 10 (30), 



provide the most stringent constraints on \k 32 \ ■ While r — > [ir\ is mediated only by A, r — > pr/ depends on the 
masses of h, H, A and H^. For sin 2 (a — /3) ~ 1 and m A ~ rn H , the branching ratio of r — > prf is suppressed because 
of the cancellation between the one-loop diagrams of A and H. Therefore, IK32I 2 is bounded most strongly by the 
r — > pr\ result for this case 3 . When tua differs from mu or when sin 2 (a — (3) is to some extent smaller than unity, 
the one-loop induced r — ► p,j process becomes important, and gives the most stringent bound on IK32I 2 of all the 
rare tau decay processes. 

The value of |K32 ax | 2 can be much larger than 100 in a wide range of the parameter region. One might think that 
such large values of |k3j| cannot be consistent with the unitarity argument for the LFV Yukawa couplings. However, 
it should be emphasized that the above figures show the contour plots for Ift^*! 2 under the rare tau decay results, 
and not for IK32I 2 . The region of | «^ 32 ax | 2 ^ 1 should be taken as the area where |«32| 2 can be as large as 0(1O~ 2 -1O~ 4 ) 
easily. It is concluded that current results of the tau LFV decays do not give any substantial upper limit on IK32I 2 
except for high tan/3 region (tan/3 > 30). 

Finally, wc comment on the bound on |/«3i| 2 , the LFV parameters for r-e mixing. Similar to T-p, mixing, we can 



discuss 



v 31 



comparing the data of r 



er?, r 



ep+p , 



eK + K and r 



e7 



listed in Table [I] with the 



3 The MSSM result approximately corresponds to this case ll4l . 
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formulas given in Ref. 24] . These formulas for r-e mixing are common with r-fi mixing except for the factor of 
|f«3i| 2 , so that difference in contours of Ik™^ 2 from those of |/«^ ax | 2 only comes from that in the data. In Table [I] 
the experimental limit for the branching ratio of r — > erj is about 1.5 times weaker than that of r — > \ir\, while that 
of t — > e7 is 5.7 times relaxed as compared to that of r — > /17. Moreover, the upper limit on Br(r~ — > e~ K + K~) is 
1.8 times stronger than that on Br(r~ — > fi~ K + K~). We have numerically confirmed that there are some regions 
where t" — > e~K + K~ can give the most stringent bound on |/t3i| 2 . Therefore, |ft3i ax | 2 is determined from one of 
t — > erj, t — > e^f and t~ — > eT K + K~ depending on parameter regions. 



4. Lepton flavor violating Higgs boson decays 

As shown in the previous section, the LFV Yukawa couplings can be tested only in the large tan/3 region by 
searching for rare tau decays. In order to cover the region unconstrained by rare tau decay results, we here consider 
LFV via the decay of the neutral Higgs bosons; i.e., <j)° — ► T^if ((f>° — h,H and A). Branching ratios for these decays 
calculated^ El El El to be 

( ™ 2 Y 

2 2/ 771 h I 1~ — T I 

Br(*-r-4-)= ^ VjT^ kl 2 r\ h ^2/ , (8) 



are 



.1 ».f si„-,„- ;i 1,, jj ••'"(' ^f") 

r(ff-»all) 

. -mi: 



Br(ff -> r^+) = ^ m ^y> |«*| a r^^V , 0) 



Br(4 = l^^gl^r r^Jau) 7 > (10) 

where r((/>° — > all) is the total width for corresponding neutral Higgs boson <p°. We here neglect terms of O{m 2 Jm 2 , ). 
Branching ratios for </>° — * t + £~ coincide with those for <fi° — > r~£f given in Eqs. JHJ, Jj|J| and (|10|l . In the following, 
we concentrate on the decays into a t-/i pair. We take the values of the SM parameters as a cm — 0.007297, 
Gp = 1.166 x 10- 5 GcV~ 2 , m z = 91.19 GeV, m T = 1.777 GeV, m M = 0.1057 GeV, m b = 4.1 GeV, m t = 174.3 GeV, 
m c = 1.15 GeV, m s = 0.120 GeV. 

A search for the LFV decays h — > T £f can give important information for extended Higgs sectors and thus for the 
structure of new physics, even when only h is found and any other direct signals for the extended Higgs sector are 
not obtained by experiments. We here evaluate the possible maximal value of the branching ratio Br(/i — > r _ /z + ) Blax 
under the results of the rare tau decay search, by inserting | «32 ax 1 2 °f Eq. J7J into the |«32| 2 in Eq. JHJ. 

In Figs. 0(a) and[2J-(b), contours of Br(/i — > r ± ^ =F ) max , which is twice of Br(h — > T~/i + ) maX) are shown in the 
tan/3-m^ plane and in the tan /3-sin 2 (a — /?) plane, respectively. The parameters are taken to be the same as those for 
Figs. ^(a) and^(b), respectively; i.e., (a) m h — 120 GeV, m H = m H± — 350 GeV and sin(a — (?) = —0.9999, and 
(b) m h = 120 GeV and m H — m A = m H± = 350 GeV, with M to be scanned from to 1000 GeV. We again show 
the excluded area from requirements of tree-level unitarity and vacuum stability as in the same way as Figs.fJ-(a) 
andEJ(b). Fot low and moderate values of tan/3 (tan/3 < 30), where rare tau decay results cannot give substantial 
upper limit on |«32 1 2 , Br(/i — ► t fi^) max can be sufficiently large. We find that the possible maximal values of the 
branching ratio can be greater than 0(1O -3 ) in a wide rage of the theoretically allowed region. 

For relatively lower tan/3 values, the experimental upper limits on IK32I 2 from rare tau decays are weaker, and 
Br(/i — > t fi^) can be sufficiently large (> 0(1O~ 3 ) for ~ 120 GeV). It is expected that a sufficient number of 
such light h can be produced at future colliders such as CERN LHC, currently planned International Linear Collider 
(ILC) and CERN CLIC. It has been pointed out that the decay process h — » r ± /i T can easily be detected at ILC 
with the luminosity of 1 ab -1 , when mh ~ 120 GeV and Br(/i — * r ± /i T ) > 0(1O~ 3 ) via the Higgsstrahlung process 
by using the recoil momentum of Z boson Therefore, the LFV search via the decay h — > t^/m* at ILC is 
complementary to that via rare tau decays at (super) B factories, and the both cover a wide region of the parameter 
space of the lepton flavor violating THDM. 

Next we discuss branching ratios for the LFV decays of heavier Higgs bosons, H/A — > T i^, using Eqs. © and 
HI U| l under the current data of LFV rare tau decays. In the THDM, there are many possible decay modes for H 
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Figure 2: Contour plots of Br(/i — *■ r ^) max , the possible maximal values for the branching ratio under the tau rare decay 
results, are shown (a) in the tan/3-mA plane and (b) in the tan/3-sin(a — 0) plane. The parameters are taken as the same as 
Figs.0-(a) and0-(b), respectively. The dark (light) shaded area indicates the excluded region by the theoretical requirements 
of vacuum stability and perturbative unitarity with f = 1 (£ = 1/2). 
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Figure 3: Contour plots of Br(_ff — ► t fj,^) iaax are shown in the tan/3-m s plane (m $ = m H — m A = m H ±) for m h = 120 
GeV and sin(a — 0) = —1 with (a) M = m 4> , (b) M — m $ /\/2, and (c) M = 0. The dark (light) shaded area indicates the 
excluded region by the theoretical requirements of vacuum stability and perturbative unitarity with £ = 1 (£ = 1/2). 



depending on the mass spectrum. In the numerical analysis, we included contributions from one-loop induced Z*y, 77 
and gg modes in addition to all the tree level modes. Those for A are one-loop induced modes of hh, hH, HH, hj, 
Hj, W^W 11 , ZZ, Zj, 77 and gg in addition to all the tree level modes. The branching ratios for H/A — » T^£f are 
sensitive to the masses of all the Higgs bosons. Here we consider the case of sin(cv — 0) = — 1 and m H = m A = rn H± 
(= m $ ). As discussed in Sec. II, the p parameter constraint is satisfied for this choice. From the 6—5-57 results, 
is taken to be greater than 350 GeV. As also discussed in Sec. II, M determines the decoupling property of 
heavier Higgs bosons. Although the branching ratios ~Br(H/A — > T^tf) are insensitive to M in the present parameter 
set, its value strongly affects the allowed parameter region under the theoretical conditions of vacuum stability and 
perturbative unitarity. Notice that couplings of H are similar to those of A for sin(a — 0) = — 1 where there are no 
HVV couplings. Hence we show the results only for the LFV decays of H below. In a general case, the branching 
ratio of H — > r ± /i T tends to be smaller than that of A — * -t^/jT due to the contribution from the modes H — > VV . 
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(a) (b) 
Figure 4: The Feynman diagram of the signal event (a), and that of the fake event (b). 



In Figs. EH-(a), El-(t>) and|31-(c), contour plots of Br(_ff — > T ± // =F ) max , the upper limit of Bi(H — > r^fji^) under the 
rare tau decay results, are shown in the tan/3-ra $ plane for M = to$, m^/y/2 and 0, respectively. As expected, the 
contours are insensitive to the values of M, and approximately the same in Figs. EJ(a),[3J(b) and|3J-(c). It is shown 
that Br(H -> r^^raax 

can be larger than 10 3 except for large tan j3 values with relatively small m $ . Therefore, it 
turns out to be no substantial upper limit on the Br(_ff — > r fJ^) in the relatively low tan (3 region (tan (3 < 20) from 
the LFV rare tau decay results. When M is smaller than m$, where the heavier Higgs boson partially receive their 
masses from the vacuum expectation value, the allowed parameter region is strongly constrained by the requirements 
of vacuum stability and perturbative unitarity. In particular, for M = (Fig. EKc)), the allowed region is limited 
only the area of around tan j3 ~ 1 and < 600 GeV. 

The extra Higgs bosons (H , A and H ± ) are expected to be searched at the LHC. The signal of gg — > H/A — > r ± /i =F 
may be detectable at LHC with high luminosity (100 fb _1 ) when Br(H/A — > T ± ^i T ) is greater than 10~ 2 for 
m H/A ~ GeV and tan/? = 45 [2^. However the rate is rapidly reduced for smaller values of tan/3 and for larger 
values of rn H , A . Further feasibility study is necessary. 



5. Search for LFV decays of Higgs bosons at a Linear Collider 

Let us consider the LFV Higgs decay h° — > r ± /i T at a LC in the situation where the heavier Higgs bosons nearly 
decouple from the gauge bosons; i.e., sin(a — 0) — — 1. The lightest Higgs boson then approximately behaves as the 
SM one. The main production modes of the lightest Higgs boson at a LC are the Higgsstrahlung e + e~ — > Z* — > 
and the W fusion e + e~ — > (VF + *j> e )(l4 y ~*i/ e ) — > h?v e v e . For a light ft with the mass ~ 120 GeV, the former 
production mechanism is dominant at low collision energies (y/s < 400-500 GeV), while the latter dominates at 
higher energies. For our purpose, the Higgsstrahlung process is useful because of its simple kinematic structure. The 
signal process is then e+e- -> Z* -> Zh° -> Zr^vF. We can detect the outgoing muon with high efficiency, and its 
momentum can be measured precisely by event- by- event. The momentum of the Z boson can be reconstructed from 
those of its leptonic £ + £~ (^ ± = e^ and /i^) or hadronic (jj) decay products. Therefore, we can identify the signal 
event without measuring r momentum directly, as long as the beam spread rate for ^/s is sufficiently low. 

Depending on the Z decay channel, the signal events are separated into two categories, jjr^fj,^ and f + ^ _ r ± /i =F . 
The energy resolution of the Z boson from hadronic jets jj is expected to be 0.3^/Ez GeV and that from £ + £~ is 
0.1y/Ez GeV. We assume that the detection efficiencies of the Z boson and the muon are 100 %, the rate of the beam 
energy spread is expected to be 0.1 % level, the muon momentum is measured with high precision and the mass of 
the lightest Higgs boson will have been determined in the 50 MeV level. We also expect that the effect of the initial 
state radiation is small for the collider energies that we consider (y/s ~ 250-300 GeV). Taking into account all these 
numbers, we expect that the tau momentum can be determined indirectly within 3 GeV for jjr fi^ and 1 GeV for 
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Let us evaluate the number of the signal event. We assume that the energy yfs is tuned depending on the mass of 
the lightest Higgs boson: i.e., we take the optimal yfs to product the lightest Higgs boson through the Higgsstrahlung 
process. (It is approximately given by yfs ~ mz + \/2m h .) The production cross section of e + e~ — > Zh° is about 
220 fb for m h = 123 GcV. Then, we obtain 2.2 x 10 5 Higgs events if the integrated luminosity is 1 ab . When | K32 1 2 
is 8.4 x 10~ 6 , about 118 events of jjr (j^ and 11 events of l + l~r ^ can be produced. 

Next, we consider the background. For the signal with the Higgs boson mass of ~ 120 GeV, the main background 
comes from e + e~ — > Zt + t~ . The number of the Zr ± /i =F event from e + e~ — > Zt + t~ is estimated about 3.6 x 10 4 . 
Although the number of the background events is huge, we can expect that a large part of them is effectively 
suppressed by using the kinematic cutsQ. The irreducible background comes from the process shown in FigQJ-(b): 
the Higgs boson decays into a tau pair, and one of the tau decays into a muon and missings (e + e~ — > Zh — > 
Zt + t~ — ► Zr ± /i T +missings). We can not distinguish the signal event h° — ► r^fi^ with the event of FigHJ-(b) when 
the muon emitted from the tau lepton carries the similar momentum to that of the parent, because it leaves the same 
track on the detector as the signal event. We refer this kind of the background as the fake signal. As examined in 
Ref. [26j , the number of the fake signal strongly depends on the precision of the tau momentum determination by the 
recoil method. We expect that it is attained with the similar precision to that of the Higgs boson mass reconstructed 
by the recoil momentum. We here take the uncertainty of the tau momentum as 3 GeV for jjr^/j^ and as 1 GeV 

for e+e-^/j^. 

Finally, we estimate the statistical significance (S/\fB) for each channel. The number of the fake events is evaluated 
in Ref. Q, which is 460 for jj^fi^ and 15 for i+fr^fi^. Therefore, when |k 32 | 2 is 8.4 x 10~ 6 with m h = 123 GeV, 
the significance can become 5.5 and 3.0 for jjr^fi^ and ^ + ^ _ r ± /i =F , respectively, taking into account the constraint 
from the rare tau decay results. The combined significance can reach to 6.3. When |«32 1 2 is 3.8 x 10~ 6 with m h = 123 
GeV, the number of the signal becomes smaller, and the combined significance amounts to be as large as 2.0. 

6. Summary 

Lepton flavor violating decays of Higgs bosons have been studied in the general THDM, in which LFV couplings 
are introduced as a deviation from Model II Yukawa interaction in the lepton sector. The model parameters are 
constrained by requirements of tree-level unitarity and vacuum stability, and also from the experimental results. The 
parameters \n3i\ 2 in LFV Yukawa interactions are bounded from above by using the current data for rare tau LFV 
decays. In the large tan/3 region (tan/3 > 30), the upper limit on \n3i\ 2 due to the rare tau decay data turns out to 
be substantial and comparable with the value predicted by assuming some fundamental theories such as SUSY. For 
smaller values of tan /3, the upper limit is rapidly relaxed, and no more substantial constraint is obtained from the 
rare tau decay results. 

We have shown that a search for the LFV decays </>° — > t if of neutral Higgs bosons ((f) = h,H and A) can 
be useful to further constrain the LFV Yukawa couplings at future collider experiments. In particular, the decays 
of the lightest Higgs boson can be one of the important probes to find the evidence for the extended Higgs sector 
even when the SM-like situation would be preferred by the data at forthcoming collider experiments. The branching 
ratio for h — > t ± // =f can be larger than O(10 -3 ) except for the high tan/3 region under the constraints from the 
current experimental data and also from the theoretical requirements. At ILC (and in case at LHC), these branching 
fractions can be tested. Therefore, we conclude that the search of LFV in the Higgs boson decay at future colliders 
can further constrain the LFV Yukawa couplings in the parameter region where rare tau decay data cannot reach. 

Note added: Recently similar work was done in Ref. [4l| on the experimental upper bound on |K3i| 2 . 
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